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ABSTRACT 

This paper discusses the out-of-plane deformation 
behavior of unsymmetric cross-ply composite plates 
compressed inplane by displacing one edge of the plate 
a known amount The plates are assumed to be initially 
flat and several boundary condition cases are consid- 
ered. Geometrically nonlinear behavior is assumed. The 
primary objectives are to study the out-of-plane behav- 
ior as a function of increasing inplane compression and 
to determine if bifurcation behavior and secondary 
budding can occur. It is shown that, depending on the 
boundary conditions, both can occur, though the charac- 
teristics are different than the pre- and post-buckling 
behavior of a companion symmetric cross-ply plate. 
Furthermore, while a symmetric moss-ply plate can 
postbuckle with either a positive or negative out-of- 
plane displacement, the unsymmetric cross-ply plates 
studied deflect out-of-plane only in one direction 
throughout tike range of inplane compression, the direc- 
tion again depending on the boundary conditions. 

miRQP Uc aQN 

A composite laminate is said to be symmetric if 
for every layer to one side of the laminate reference sur- 
face with a specific thickness, specific material proper- 
ties, and specific fiber orientation, there is another layer 
the identical distance on the opposite side of the refer- 
ence surface with identical thickness, material proper- 
ties, and fiber orientation. If the laminate is not 
symmetric, then it is referred to as an unsymmetric lam- 
inate [1]. The terminology nonsymmetric and asymmet- 
ric are also used. There are special cases of unsymmetric 

laminates, v/uc sulu Xo ou qnm j u t itiw luc cross ply lami- 
nate, whereby the laminate consists of an even number 
of layers with the principal material directions of the 
layers alternating at 0° and 90° to the laminate axes. An 
antisymmetric angle-ply laminate has layers oriented at 
+0 to the laminate coordinate axes on one side of the 
reference surface and layers oriented at -0 to the lami- 
nate axes on the other side of the reference surface [2]. 
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If a laminate is unsymmetric, then the response charac- 
teristics are somewhat unconventional due to the cou- 
pling of inplane forces with bending deformations and 
holding moments with inplane deformations, couplings 
not present in symmetric laminates. The level of asym- 
metry, and therefore the magnitude of the resulting elas- 
tic couplings, is determined by the lamination sequence. 
The higher the level of asymmetry, the higher the mag- 
nitude of the elastically-coupled responses. The level of 
asymmetry of a particular laminate may be evaluated by 
considering the laminate stiffness matrix, namely, the 
ABD matrix. Specifically, the B matrix, a submatrix of 
the ABD matrix, determines the level of coupling 
between the inplane and out-of-plane responses, also 
known as bending-stretching coupling. Laminates may 
be tailored to exhibit certain elastic couplings to take 
advantage of this unique capability. Unfortunately, tradi- 
tional elevated-temperature cure techniques make it dif- 
ficult to produce an unsymmetric laminate that does not 
change shape as it is cooled from the cure temperature 
[3]. As a result, flat unsymmetric laminates axe difficult 
to manufacture. One approach is to form an unsymmet- 
ric laminate by bonding together two symmetric lami- 
nates using a room-temperature adhesive. However, due 
to the limitations of room-temperature adhesives and the 
added cost of manufacturing, this approach is of limited 
applicability. An alternative manufacturing method that 
may hold promise for the future is the electron beam 
curing process. Due to the low temperatures associated 
with tius process, thermally-induced shape changes may 
be minimized or eliminated altogether. An even more 
intriging idea would be to use nanotechnology to tailor 
the properties through the thickness of a flat plate to be 
unsymmetric. 

Once flat, or curved, unsymmetric laminates 
can be manufactured with a degree of dimensional fidel- 
ity, attention will then focus on understanding the 
response of unsymmetric laminates to mechanical, and 
perhaps, thermal loads. The unique ability to couple 
inplane response with out-of-plane response will make 
unsymmetric laminates candidates for advanced struc- 
tural tailoring concepts and possible solutions to config- 
uration changes needed for adaptive structures not 
achievable by other means. Before exploring those 
applications, however, it will be necessary to fully char- 
acterize the response of unsymmetric laminates to sim- 
ple loadings, categorize their behavior, and understand 
their limits. This paper is directed at these goals and will 
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consider the deformation response of flat unsymmetric 
laminates to a simple inplane compressive loading. This 
loading is considered because it is one of the simplest 
mechanical loads to model and one that can be consid- 
ered in a laboratory setting. For symmetric laminates 
such a load would be associated with the prebuckling, 
buckling (bifurcation), and postbudding behavior. For 
unsymmetric laminates, which, due to the elastic cou- 
plings, can begin to deform out of plane as soon as they 
are compressed inplane, it is not clear if the concepts of 
prebuckling, budding, and postb tickling are valid. That 
very point is part of this investigation. Leissa [4] argued 
that a flat unsymmetric laminate will remain flat under a 
compressive inplane loading if out-of-plane rotations 
are restricted an all four edges, i.e., the plate is clamped. 
This occurs because the clamped conditions provide the 
moments necessary to counter the moments induced by 
the inplane compression which make the plate deform 
out of plane. If this flat state exists, then classic preb- 
uckling, buckling, and postb uckling responses are a pos- 
sibility. Leissa [5] went on to study the geometrically 
linear responses of four types of unsymmetric laminates 
subjected to uniform and linearly-varing inplane loads, 
including uniaxial, biaxial, and shear loads. Twelve sets 
of boundary conditions were considered. The result of 
the study was a tabulation of the center deflections of 
the many cases, including the case of no deflection for 
the conditions discussed in ref. 4. The present paper 
could be considered an extension to ref. 5 in that the 
compressive loading discussed here is one of the load- 
ings Leissa [5] considered, and plate response into the 
geometrically nonlinear range is considered to further 
address the issues of prebuckling, buckling, and post- 
buckling. However, only cross-ply laminates are consid- 
ered here and fewer boundary conditions are 
investigated, as the response of the laminates becomes 
complex when in the postbuckled state. 

The following section describes the specific 
problem studied. The loading, boundary conditions, 
coordinate system, geometry, and nomenclature are 
described As the numerical results are determined using 
the commercial code ABAQUS [6], finite-element con- 
siderations are also discussed. The section following 
that discusses numerical results. 

PROBLEM DESCRIPTION AND ANALYSIS 

DETAILS 

The geometry and nomenclature of the prob- 
lem studied are shown in fig. 1. Rectangluar plates of 
dimensions L by W are considered, though one square 
plate will be discussed. The long dimensions of the plate 
are referred to as the sides and the short dimensions the 
ends. The x-coordinate originates at one end of the plate, 
and the y-coordinate is centered between the two sides. 


The standard principal material coordinate system for 
the layers is located within this coordinate system, as 
shown. The fiber, or 1, direction makes an angle 0 rela- 
tive to the +x-axis. The displacements in the x- f y-, and 
^-coordinate directions are denoted by u, v, and w, 
respectively. The rotations about the coordinate axes are 
denoted as and and fig. 1 shows the sense of 
these rotations as defined in ABAQUS. The compres- 
sive loading is applied by displacing the right end of the 
plate a known amount U to shorten die plate. The asso- 
ciated load is denoted as P. Such a loading can be pro- 
duced in the laboratory, and since it makes sense to 
begin any experimental investigation of the deformation 
response of unsymmetric laminates by considering 
small-scale laboratory size specimens, the plates consid- 
ered are 0.508 m long by 0.406 m wide (20.0 in. by 16.0 
in.). The plates are asumed to be constructed of eight- 
layers of graphite-epoxy and measure L016 mm (0.04 
in.) thick. The properties of a layer are assume to be 

E 1 ~130.0GPa (18.85 Msi) 

E 2 *=9.70GPa (1.407 Msi) 

G 12 =5.00 GPa (0.725 Msi) 
v i 2 *0300 

h=«.1270mm(0.005in) 

Out-of-plane displacement vs. endshortening results 
will be presented at the three locations labeled /, //, and 
III in fig. 1. Locations / and III are quarter locations, 
while location II is at the center. Plots of the deformed 
geometry will be presented for specific eadshartening 
levels. Additionally, the relationship between endshort- 
eaing and the associated load P will be presented. This 
relation is often used for problems such as discussed 
hoe to describe the inplane stiffness of the plate as the 
response changes from prebuckling to postbudding. 
Three combinations of boundary conditions will be pre- 
scribed along the ends and sides, and the result* will be 
labeled SS-SS, CL-SS, and CL-CL to designate the par- 
ticular boundary condition combination being consid- 
ered. The terminology SS is contracted notation for 
simply-supported conditions and CL is a contraction for 
clamped conditions. The contraction to the left of the 
hyphen designates the prescribed boundary condition 
along the ends, while the contraction to the right of the 
hyphen designates the prescribed boundary condition 
along tiie sides. What is important to realize from the 
onset of the discussions is that the sides of the plates are 
assumed to be constrained so that iH) there. This might 
be considered a somewhat arbitrary and unrealistic con- 
dition, but it is also unrealistic to consider plates that 
might be joined with other elements to form a complete 
structure to have no restrictions on the v-displacements 
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along the sides. Rather than look at a range of elastic 
restraints, adding another variable to the problem, the 
case of full restraint is considered. This issue is impor- 
tant because, through Poisson effects, inplane compres- 
sion in the jr-direction can cause the plate to expand in 
the ^-direction. With v=0 along the sides, Poisson 
expansion is restricted by an inplane compression force 
in the ^-direction. With the bending-stretching coupling 
present in an unsymmetric laminate, this inplane farce, 
along with the applied inplane force, arc factors. 

Ninc-noded, two-dimensional, geomctrically- 
nonlinear, shear-deformable SSRS shell elements are 
used to model the plates. The plates are compressed stat- 
ically by increasing the endshortening using very small 
increments of U. At each value of U the static stability 
of the plate is checked by examining the eigenvalues of 
the tangent stiffness matrix. A zero or negative eigen- 
value signals a statically unstable laminate configura- 
tion. At the value of U where the static solution ceases 
to be stable, some iteration is then done to adjust U so 
that the value of (/just at the onset of the instablity can 
be closely determined. For a value of U slightly greater 
than this onset value, and for this value of U held fixed, 
a dynamic analysis, using a small amount of damping, is 
then initiated. For a statically unstable condition, simply 
initiating a dynamic analysis should be enough to allow 
the plate to begin to move toward a stable equilibrium 
configuration. However, for some unstable conditions, 
this requires integrating forward in time for long periods 
of time, as the plate may move slowly from its unstable 
configuration. Accordingly, as an alternative to simply 
initiating the dynamic analysis, the plate is given a very 
small pressure pulse perpendicular to the plate at time 
zero of the dynamic analysis. The dynamic analysis is 
carried out for small pressure pulses in both the +z and - 
z directions (see fig. 1). In general, if pressure pulses are 
used, the plate moves to the same configuration inde- 
pendent of the sign of the pressure pulse. As the plate 
motion is decaying around this unique configuration, the 
decaying being due to the light damping, a new static 
analysis is inititated using a decayed dynamic configu- 
ration as a starting point The static analysis then con- 
verges within a reasonable number of iterations to what 
is interpreted to be the statically stable configuration for 
that value of enshortening. This numerical scenario cor- 
responds to a displacement-controlled compression test 
in the laboratory. Ten elements in each direction, for a 
total of 100 elements, are used to model the plates. To 
test the robustness of this mesh, for several specific test 
cases the mesh was doubled in each direction, for a total 
of 400 elements. All aspects of plate response were 
identical for the two meshes, including the dynamic sta- 
bility analysis. Of course, the dynamic stability analysis 


took considerably more time with 400 elements than 
with 100 elements, so the 100 element mesh was used. 

NUM ERICAL RESULTS 

Although the goal of this paper is to investigate 
the response of unsymmetric composite plates, results 
will first be presented for an initially fiat CL-SS [0^ 
902 ]g symmetric cross-ply laminate. Because it is sym- 
metrically laminated and flat, this plate remains flat 
when initially loaded, and then exhibits classic bifurca- 
tion behavior. The response for this situation is well 
known, so relatively little new information will be pre- 
sented regarding this particular plate. Rather, the results 
presented serve as a baseline for comparing the behavior 
of the initially flat unsymmetric cross-ply lam- 

inate to be discussed. Also, by considering the symmet- 
ric cross-ply plate, the dynamic stability analysis 
implemented in ABAQUS is applied to a problem that is 
somewhat familiar, and die characteristics of the 
approach can be studied. The approach can then be 
applied to the unsymmetric laminate with some experi- 
ence and confidence. 

[ 0 2 / 902 ] s Symmetric Cross-Ply Laminate 

Figure 2 presents the normalized load vs. end- 
shortening relation, the normalized out-of-plane dis- 
placement vs. endshortening relations, and deformed 
plate configurations at selected values of endshortening 
for the CL-SS [0 2 /902]s symmetric cross-ply laminate as 
the endshortening increases from zero to a normalized 
value of 10. The normalizing factor Ucr is the critical, 
or buckling, value of endshortening computed by 
assuming a linear, non-shear-deformable eigenvalue 
solution to the classic SS-SS uniaxial buckling problem 
[2]. The term ‘classic SS-SS’ implies that the bound- 
aries are flee to move inplane (v*0) and are only con- 
strained from o'lt-of-nlane deformation (w=*0). The 
value of Ucr for the specific plate here is 0.00516 mm 
(0.000302 in.) and this normalizing factor was chosen 
simply because it is easily calculated. The CL-SS 
boundary conditions increases the critical aidshortening 
with respect to the classic value, and therefore, a nor- 
malized aidshortening value of 10 does not imply that 
the plate being discussed is compressed uniaxially to 10 
times its own critical value. Strain levels in the plate 
were checked at U/Ucr=10 and they were not particu- 
larly high. The normalizing factor Per used to normalize 
P is the load associated with Ucr. 

The three figures for the out-of-plane displace- 
ments at points I, II, and HI show that as the endshorten- 
ing is increased from zero, point A, the plate remains 
flat At endshortening values past point B y the plate 
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remains flat and is in equilibrium but it is unstable, as is 
well known. At a point just past point B in the flat con- 
figuration, the dynamic analysis is initiated. Since the 
plate is very close to being in a stable condition, simply 
integrating the equations of motion with no small pres- 
sure pulse results in a large number of time integration 
steps far the plate to converge to a stable equilibrium 
configuration. A small pressure pulse, e.g., 689 Pa (0.1 
psi), applied in a ramp-up/ramp-down fashion, hastens 
the motion to a stable configuration characterized by a 
small out-of-plane displacement, as indicated by point C 
in the three figures and as shown by the leftmost 
deformed configuration at the bottom of fig. 2. This 
deformed configuration has one halfwave in the loading 
direction and one halfwave perperdicular to the loading 
direction. This, of course, is very much like the first lin- 
ear buckling mode for this plate, but at point C the plate 
is in a postbuckled state. If the pressure pulse is positive, 
the plate moves to configuration C, and if the pulse is 
negative, the plate moves to configuration C\ Increas- 
ing the endshortening toward point A the single half- 
wave configuration deepens, as shown at the bottom of 
the figure. Past point D equilibrium configurations exist 
but they are again unstable. Initiating a dynamic analy- 
sis just past point D results in the plate moving to con- 
figuration E. The plate could also move to configuration 
E\ If after point B the plate response is given by path 
C’D\ then initiating a dynamic analysis at D’ would 
result in the plate also moving to configuration E or£.’ 
As can be seen, the configurations associated with E and 
E* are represented by two halfwaves in the loading 
direction and one halfwave perpendicular to the loading 
direction. This change of configuration is often referred 
to as secondary buckling and is discussed for composite 
plates by Tiwari and Hyer [7]. Increasing the endshort- 
ening to point F deepens the configuration of point E. 
After point 2?, the slope of the normalized load vs. rad- 
shortening relation, in the upper left of fig. 2, decreases, 
and the relation becomes slightly nonlinear, sullening as 
configuration A D\ is approached. The change to 
configuration E or E* is accompanied by a slight drop in 
load. Beyond point E t and the slope of the relation 
decreases again. 

[<V904]t Unsymmetric Cross-Ply Laminates 

If the eight layers of the [O2/9OJS symmetric 
cross-ply laminate are rearranged to form a [(V90JT 
unsymmetric cross-ply laminate and the numerical anal- 
ysis of the previous section far CL-SS boundary condi- 
tions is repeated, the normalized load vs. endshortening 
relation, die normalized out-of-plane displacement vs. 
endshortening relations, and deformed plate configura- 


tions at selected values of endshortening shown in fig. 3 
result The format of fig. 3 is identical to fig. 2 to aid in a 
direct comparison of plate responses. A quids, glance at 
the out-of-plane displacements at locations /, A and HI 
shows that unlike the symmetric cross-ply 

case with multiple stable solution paths (BCD, B'C'D\ 
EF and E'F\ there is only one stable solution path. As 
the endshortening is increased from zero, the plate 
deforms in a stable fashion out of plane in the negative 
z-diiection a very small amount Standard experimental 
techniques probably could not detect any displacement, 
but the displacements are nonzero and the configuration 
at B is one with two halfwaves in each direction, as 
shown at the lower left in fig. 3. Far endshortening val- 
ues past point B y the plate is unstable and the dynamic 
analysis results in the plate moving to the configuration 
given by point C. Point C represents a large increase in 
the out-of-plane displacement relative to point B y but the 
configuration at Point C, shown at the bottom of the fig- 
ure, has a configuration like point B . Many attempts 
were made to determine if the gap between points B and 
C is actually closed. This does not seem to be the case. 
With the finite-element analysis it is not possible to have 
the endshortening values increase continuously. How- 
ever, even using fine increments in endshortening 
beyond point B y there appears to be a jump from 21 to C. 
This suggests that this particular unsymmetric laminate, 
with these particular boundary conditions, does not 
exhibit bifurcation. Increased endshortening to point D 
further deepens the halfwave configuration. For an end- 
shortening value just past point A a dynamic analysis is 
again initiated and stable configuration E results, inde- 
pendent of the sign of the small pressure pulse. The con- 
fig uation at A shown at the bottom of the figure, is not 
symmetric about the x=U2 line, as can be seen by com- 
paring the out-of-plane displacement at location I with 
that at location III. Increased endshortening beyond con- 
figuration E y which exhibits negative out-of-plane dis- 
placement at all three locations, though they are 
relatively small at location III , results in a gradual tran- 
sition to the configuration at point F. This configuration, 
as shown at the lower right, has a large negative out-of- 
plane displacement at location /, a smaller negative out- 
of-plane displacement at location //, and positive out-of- 
plane displacement at location III. The configuration has 
two halfwaves in the loading direction, but one that 
skewed to one end of the plate rather than being anti- 
symmetric with respect the the x=L!2 location. 

If the boundary conditions on the ends of the 
plate are changed to be simply supported, then the 
responses shown in fig. 4 result As can be seen, chang- 
ing of the boundary conditions has a profound influence 
on die response of the plate. Specifically, as die end- 
shortening is increased from zero the plate deforms out 
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of plane a significant amount This is in contrast to the 
CL-SS case previously discussed in fig. 3, where the 
plate remained relatively flat for low end shortening lev- 
els. For the SS-SS case, as seen in the out-of-plane dis- 
placement vs. endshortening relations at locations 7, //, 
and III in fig. 4, up to point B the plate deforms out of 
plane in a single halfwave configuration. At point 2?, the 
load vs. endshortening relation decreases slope slightly. 
As the endshortening is increased, the plate configura- 
tion begins to change. When the aidshortening reaches 
the level of point C, the plate appears to have two half- 
waves in the loading direction. A close look at the out- 
of-plane displacements at location 77, however, reveals 
that the displacement is not zero in the center of the 
plate. At point C the plate becomes unstable, and a 
dynamic analysis is used to affect the transition to point 
D. At point D the deformed configuration looks like that 
at point C, but the out-of-plane dispacement at location 
77 is positive, while the magnitude of the out-of-plane 
displacements at locations 7 and in has decreased rela- 
tive to point C. A further increase in endshortening to 
point E does not produce any configuration changes. It 
is clear there is nothing that resembles bifurcation 
behavior near point B in fig. 4, rather the plate configu- 
ration makes a smooth continuous shape change as the 
endshortening increases from A to B to C. Furthermore, 
the plate deflects out of plane in a direction opposite to 
the CL-SS case in fig. 3. 

If instead of enforcing simply-support bound- 
ary conditions on all four sides, clamped conditions are 
enforced, the responses illustrated in fig. 5 occur. As the 
endshortening is increased from zero, the plate appears 
to remain flat According to Leissa [4], this should be 
the case, though the finite-element results indicate there 
are very small out-of-plane displacements. As they are 
many orders of magnitude less than the out-of-plane dis- 
placements along path AB in fig. 3, these are felt to be a 
result of numerical round-off in the finite-element anal- 
ysis. Increasing the endshortening past point B results in 
a statically unstable condition. The dynamic analysis 
leads to a statically stable equilibrium configuration 
whereby the plate is deformed out of plane in what 
appears to be a two halfwave configuration, as shown by 
the inset of the deformed geometry of point C. However, 
the configuration is not exactly a two halfwave one, 
rather there is a small out-of-plane displacement at the 
center of the plate, location 77, as seen in the upper por- 
tion of fig. 5. With increasing endshortening past point 
B, the displacement at location 77 is at first negative, 
then positive, though it remains small. At locations /and 
III the displacments are much larger and have opposite 
signs. This psuedo-halfwave configuration remains 
throughout the range of endshortening studied. There is 
no secondary buckling behavior, as was the case for the 


CL-SS and SS-SS boundary conditions of figs. 3 and 4. 
As was the case earlier, changing the boundary condi- 
tions relative to the CL-SS case of fig. 3, in this case 
making the sides clamped, has an important effect on 
plate response. In particular, it would appear that the 
plate exhibits bifurcation behavior at 2?, much like a flat 
symmetrically laminated plate of fig. 2. However, 
instead of bifurcating to positive and negative out-of- 
plane displacement branches, only one out-of-plane dis- 
placement branch occurs. 

The psuedo-halfwave response throughout the 
endshortening range for the CL-CL boundary condition 
shown in fig. 5 is a result of the rectangular geometry of 
the plate. If the sidelength of the plate in fig. 5 is 
decreased so the plate is square, the response shown in 
fig. 6 results. As the endshortening is increase from zero 
to point 2?, the plate remains flat, to within the numerical 
accuracy of the analysis. At point B the flat configura- 
tion becomes unstable. The dynamic analysis is initiated 
for a fiat configuration just beyond point B and the con- 
figuration of point C results. This configuration has a 
single halfwave in both directions. This configuration 
remains as the endshortening is increased to point D. At 
point D tire single halfwave configuration becomes 
unstable and the dynamic analysis is again initiated. The 
plate assumes a configuration with two halfwaves in the 
loading direction and one halfwave perpendicular to that 
direction. Noting the nonzero displacement at location 
77, it is seen that the configuration is not truly a two half- 
wave one. However, as the endshortening is increased 
beyond point £, the out-of-plane displacement of the 
center of the plate approaches zero and the displace- 
ments at locations 7 and III approach being equal and 
opposite. For this specific configuration, bifurcation 
occurs at point B , and secondary buckling occurs at 
point D. The bifurcation results in just one out-of-plane 
solution branch. 

SvwiARY 

This study has shown that the behavior of unsymmetri- 
cally laminated cross-ply plates subjected to inplane 
compression can exhibit a range of out-of-plane defor- 
mation behavior, depending on the level of compres- 
sion, boundary conditions, and geometry. For this 
particular study the condition v=0 was enforced along 
the sides, further contributing to the response. Except 
for the CL-CL boundary condition, the initially flat 
plates did not remain flat of low levels of compressive 
loading, though for the CL-SS case the out-of-plane dis- 
placements were very small and could be interpreted as 
zero in an experiment For the CL-SS case the very 
small out-of-plane displacement configuration did 
become unstable and an alternative stable configuration 
was found which did have large out-of-plane displace- 
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merits. Secondary buckling could occur, resulting in 
another stable equilibrium configuration. 

ArKNOWIFPfiEMENTS 
The work discussed herein was supported by grant 
NAG-1-2298 from the Mechanics and Durability 
Branch of the NASA-Langley Research Center. Dr. 
Damodar Ambur was the grant monitor. 

REFERENCES 

1- Hyer, M.W. Stress Analysis of Fiber-Reinforced Com- 
posite Materials , WBC/McGraw-Hill, New Yoric, 1998, 
p. 296 

2 • Jones, R.M., Mechanics of Composite Materials, 2nd 
Ed., Taylor-Francis Publishers, Philadelphia, 1999 

3 - Dano, M.-L. and M.W. Hyer, “Thermally-Induced 
Deformation Behavior of Unsymmetric Laminates,” Int 
J Solids Struct, 35 (17), pp. 2101-2120, 1997 


4 - Leissa, A.W., “Conditions for Laminated Plates to 
Remain Flat under biplane Loading,” Composite Struc- 
tures, 6, pp. 261-70, 1986 

5 - Qatu, M.S. and A.W. Leissa, “Buckling or Trans- 
verse Deflections of Unsymmetrically La min a t ed Plates 
Subjected to Inplane Loads,” AHA J., 31, 1, pp. 189-94, 
1993 

6 - ABAQUS/STANDARD. Theory Manual. Hibbitt, 
Karlsson, & Sorensen, Inc., Pawtucket, Rhode Island, 
1998 

7 - Tiwari, N. and M.W. Hyer, “Secondary Buckling of 
i aminatftH Composite Plates,” AIAA Paper 2000-1457, 
presented at 41st AIAA/ASME/ASCE/AHS/ASC Struc- 
tures, Structural Dynamics, and Materials Conference, 
Atlanta GA, April 2000, also AIAA J, to be published 


y, V.-<|>y 



Fig. 1 - Problem description, geometry, end nomenclature 
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Fig. 5 - Response Characteristics of CL-CL [04/90 4 ] x Plate 
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Fig. 6 - Response Characteristics of Square CL-CL [04/90 4 ] T Plate 
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